Isotropic three-dimensional gap in the iron-arsenide superconductor LiFeAs from 

directional heat transport measurements 



M. A. Tanatar, 1 '^ J.-Ph. Reid, 2 S. Rene de Cotret, 2 N. Doiron-Leyraud, 2 F. Laliberte, 2 E. Hassinger, 2 J. Chang, 2 
H. Kim, 1,3 K. Cho, 1 Yoo Jang Song, 4 Yong Seung Kwon, 4 R. Prozorov, 1,3 and Louis Taillefer 2, 5, [j] 

1 The Ames Laboratory, Ames, Iowa 50011, USA 
Departement de physique & RQMP, Universite de Sherbrooke, Sherbrooke, Canada 
J Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011, USA 
J * Department of Physics, Sungkyunkwan University, 
Suwon, Gyeonggi-Do 440-74-6, Republic of Korea 
Canadian Institute for Advanced Research, Toronto, Ontario, Canada 
(Dated: January 19, 2013) 

The thermal conductivity k of the iron-arsenide superconductor LiFeAs (T c ~ 18 K) was mea- 
sured in single crystals at temperatures down to T ~ 50 mK and in magnetic fields up to H = 17 T, 
very close to the upper critical field H C 2 — 18 T. For both directions of the heat current, parallel 
and perpendicular to the tetragonal c axis, a negligible residual linear term n/T is found as T — > 0, 
revealing that there are no zero-energy quasiparticles in the superconducting state. The increase 
in k with magnetic field is the same for both current directions and it follows closely the depen- 
dence expected for an isotropic superconducting gap. There is no evidence of multi-band character, 
whereby the gap would be different on different Fermi-surface sheets. These findings show that 
the superconducting gap in LiFeAs is isotropic in 3D, without nodes or deep minima anywhere on 
the Fermi surface. Comparison with other iron-pnictide superconductors suggests that a nodeless 
isotropic gap is a common feature at optimal doping (maximal T c ). 



PACS numbers: 74.25.fc, 74.20.Rp,74.70.Xa 

Because the structure of the superconducting gap as 
a function of direction reflects the pairing interaction, it 
can shed light on the nature of the pairing mechanism. In 
the iron pnictides, the experimental situation in this re- 
spect remains unclear and so far suggests the lack of any 
universal picture. Several studies agree on the existence 
of nodes in the superconducting gap of the low-T c materi- 
als KFe 2 As 2 [HE] and LaFePO [3H2]. In BaFe 2 As 2 -based 
superconductors, signatures of nodal behavior were ob- 
served in heavily K-doped samples [5] and in P-doped 
compounds [7], while in Co- and Ni-doped compounds 
the superconducting gap shows nodes only away from 
optimal doping (maximal T c ) [8HIU]. 

The material LiFeAs may prove important in the study 
of iron-based superconductivity because it is stoichiomet- 
ric, and so can in principle be made with low levels of 
disorder, and it has a relatively high T c . The Fermi sur- 
face of this material has four (or five) sheets, two electron 
pockets centered near the M-point of the Brillouin zone, 
and two (three) hole pockets centered around the T-point 
[llj . ARPES measurements for k z =0 found an isotropic 
in-plane superconducting gap whose magnitude on the 
electron sheets, A e , is approximately two times larger 
than on the hole sheets, A/, [12] . Specific heat [13], pen- 
etration depth [HI [TS] and lower critical field [TB] mea- 
surements were interpreted in terms of a fully isotropic, 
k- independent gap A(fc), with A e ~ 2A^. However, none 
of these studies has directional resolution to locate out of 
plane nodes, such as three dimensional (3D) nodes found 
in the under- and over-doped Co-Bal22 [10l [J7] 

In this Letter, we report a study of the 3D supercon- 



ducting gap structure of LiFeAs using thermal conduc- 
tivity, a bulk probe used previously to locate directions 
of gap nodes in heavy-fermion [18] and iron-pnictide [TO] 
superconductors. We found that for directions of heat 
flow parallel and perpendicular to the tetragonal c axis 
the thermal conductivity of LiFeAs closely follows expec- 
tations for a single isotropic superconducting gap, with 
no evidence of nodes or deep minima in any direction on 
any part of the Fermi surface. 

Experimental- Single crystals of LiFeAs were grown 
in a sealed tungsten crucible using a Bridgeman method 
|19j . and stored in sealed ampoules. Immediately af- 
ter opening the ampoules, samples for in-plane resis- 
tivity, Seebeck and thermal conductivity measurements 
were cleaved and shaped into parallel bars (1 — 2) x 
(0.3 - 0.5) x (0.05 - 0.1) mm 3 (a x b x c). Silver 
wires were soldered to the samples [20] . yielding low- 
resistance contacts (~ 100 /ifl). Samples for inter-plane 
resistivity and thermal conductivity, with dimensions 
(0.5 — 1) x (0.5—1) x (0.1 — 0.3) mm 3 , were measured using 
a two-probe technique [TU1 [Uj ■ , with contacts covering 
the whole a6-plane area of the sample. After contacts 
were made, the samples were covered with Apiezon N 
grease to prevent oxidation. 

The thermal conductivity k was measured in a stan- 
dard one-heater-two thermometer technique. In both in- 
plane (k q ) and inter-plane (k c ) heat transport measure- 
ments, the magnetic field H was applied along the [001] 
tetragonal c axis. Measurements were done on warming 
after cooling in a constant field from above T c , to ensure 



2 





3.0 r 




2.5 




2.0 




1.5 




1.0 


E 




o 




CM 


0.5 








0.0 


E 




25 r 


\- 






X 


0.20 




0.15 




0.10 




0.05 




0.00 




a axis 


4T 




1 T: 




0.05 T 








0.000 



0.00 0.01 0.02 0.03 0.04 0.00 0.01 0.02 0.03 0.04 



T 2 (K 2 ) 



T 2 (K 2 ) 



FIG. 1: Thermal conductivity of LiFeAs as a function of 
temperature, plotted as k/T vs T 2 , for a heat current in the 
basal plane (left panels) and along the tetragonal c axis (right 
panels), measured for different values of the magnetic field H 
as indicated. Solid lines are linear fits used to extrapolate the 
residual linear term kq/T at T = 0, plotted vs H in Fig. 2. 



a homogeneous field distribution. 

Temperature dependence. - The thermal conductivity of 
LiFeAs is displayed in Fig. [T] for different magnetic fields 
up to 17 T. The linear fits show that the data below 0.2 
K are well described by the function k/T = a + bT 2 . The 
first term, a = Ko/T, is the residual linear term, entirely 
due to electronic excitations [35] . The second term is due 
to phonons, which at low temperature are scattered by 
the sample boundaries. 

The magnitude of the residual linear term is extremely 
small. For both directions of heat flow, n /T ~ 5 /iW 
/ K 2 cm. These values are within the absolute accuracy 
of our measurements, approximately ± 5 /jW / K 2 cm 
[23), [25] . Therefore, our LiFeAs samples exhibit a negligi- 
ble residual linear term for both in-plane and inter-plane 
directions. Comparison with the normal-state conduc- 
tivity «n/T, estimated using the Wiedemann- Franz law 
- Kn/T — Lq/pq where Lq = (7r 2 /3)(fcB/e) 2 - applied 
to the extrapolated residual resistivity p (see Fig. 3), as 
discussed in Ref. ITO] gives a ratio (kq/T) / (k^ /T) ~ 1 % 
(0.1 %) for flow parallel (perpendicular) to the c axis. 

These kq/T values are much smaller than theoreti- 
cal expectation for a nodal superconductor (for a gap 
without nodes, kq/T = [22]). For a quasi-2D rf-wave 
gap, with four line nodes along the c axis, the resid- 
ual linear term is given, in the clean limit, by kq/T = 
(/cg/6c)(fci7i>i?/Ao), where c is the interlayer separation, 
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FIG. 2: Residual linear term ko/T in the thermal conductiv- 
ity of LiFeAs as a function of magnetic field H (applied along 
the tetragonal c axis) , plotted on scales normalized to the nor- 
mal state. Kr-j/r is the normal-state conductivity estimated 
from the Wiedemann- Franz law (see text); H C 2 is the upper 
critical field in the T — limit (see Fig. 3). The same field de- 
pendence is observed for the two directions of heat flow, along 
(J || c) and perpendicular (J || a) to the c axis. This isotropic 
behaviour is very similar to that of standard isotropic s-wave 
superconductors, as in the clean Nb and the dirty InBi shown 
here (reproduced from Ref. I35[) . For comparison, we also re- 
produce data for the ci-wave (nodal) superconductor Tl-2201 
[36] and the multi-band s-wave superconductor NbSe2 [25] . 



Uf and vf the Fermi wavevector and velocity at the node, 
respectively, and A the gap maximum [22, 26--28J. Tak- 
ing c = 6.36 A, v F = 1 eV A = 1.5 x 10 5 m/s QUCIH], 
and a typical Fermi wavevector for electron sheets of the 
Fermi surface, k F = 0.2 (ir/a) = 0.16 A" 1 [30], we get 
Ko/T ~ 140 /iW / K 2 cm, assuming a weak-coupling 
Aq = 2.14 kftTc, not far from the experimentally deter- 
mined gap 31 . This is at least 20 times larger than the 
value extracted from our fits to the k/T vs T data. In 
those materials where universal heat transport has been 
verified, proving the presence of a line node in the gap, 
the measured value of kq /T is in good quantitative agree- 
ment with theoretical expectation [15 ] I28 | 1521 133] . Thus 
we can safely conclude that the gap in LiFeAs does not 
contain a line of nodes anywhere on the Fermi surface. 
Importantly, the fact that kq/T ~ for both K a and k c 
rules out not only vertical but also horizontal line nodes, 
including those away from the k z = plane. 

Field dependence. - Our zero- field data show that there 
are no zero-energy quasiparticle excitations in LiFeAs, 
and therefore no nodes in the gap structure anywhere 
on the Fermi surface. By applying a magnetic field, we 
can now investigate quasiparticles at energies above zero. 
In a type-II s-wave superconductor, a field applied per- 
pendicular to the heat flow promotes heat transport by 
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allowing tunneling between the quasiparticle states local- 
ized in the core of adjacent vortices 24 . The stronger 
the field, the closer the vortices, exponentially favouring 
the tunneling process, controlled by the ratio of coher- 
ence length £o to inter- vortex separation [231 ■ For a 
full isotropic gap, this yields an exponential growth in 
K vs H, as shown in Fig. 2 for Nb. Now if the gap is 
depressed on some region of the Fermi surface - either 
by being smaller on one sheet (multi-band character) or 
by having a strong angle dependence leading to a deep 
minimum in some k direction (gap anisotropy) - the tun- 
neling will be enhanced, since £o cx vf/Aq will be longer. 
This in turn will enhance the thermal conductivity at low 
field, as observed for example in the multi-band s-wave 
superconductor NbSe 2 [35] (see Fig. 2), or in the highly 
anisotropic s-wave superconductor LuNi 2 B 2 C [33]. 

In Fig. 2, we show the field dependence of kq/T in 
LiFeAs, obtained by extrapolating the in-field n/T vs T 
data of Fig. 1. Both axes of the plot are normalized to the 
respective normal-state value, kq /T is measured relative 
to the normal-state residual conductivity k^/T = L /p , 
with the residual resistivity po obtained by extrapolat- 
ing p(T) to T = (see Fig. 3). Note that in the ratio 
(kq/T) /(kn /T) the usual uncertainties in the geometric 
factors of the samples cancel out, since heat and charge 
transport are measured using the same contacts. The 
only uncertainty lies in the T = extrapolation of n/T 
to get k /T (well below ± 10 %; see Fig. 1) and of p(T) 
to get po (of order ± 20-30 %; see Fig. 3). The field 
axis in Fig. 2 is measured relative to the T — upper 
critical field H c2 (0) ~ 18 T, obtained by smoothly ex- 
trapolating H vs T c data to T c = 0, where T c is detected 
in thermopower measurements on LiFeAs (see Fig. 3). 
The value H C 2(0) — 18 T is consistent with tunnel-diode- 
resonator measurements on the same batch of crystals 
[37]. 

In Fig. 2, the field dependence of k /T in LiFeAs is seen 
to be isotropic, slow at low H and rapid as H approaches 
H C 2- This upward curvature of kq/T vs H is typical 
of isotropic s-wave superconductors like Nb (clean limit) 
and InBi (dirty limit), as shown in Fig. 2. It is opposite to 
the field dependence expected for a gap with nodes [2"2"] . 
as illustrated in Fig. 2 with data for the <i-wave cuprate 
superconductor Tl-2201 [36 . In this case, the Doppler 
shift of delocalized quasiparticle excitations (not confined 
to the vortex cores) yields a rapid initial rise [38]. 

The in-field data not only confirms the absence of 
nodes in the gap of LiFeAs, it also shows that the gap is 
isotropic in 3D, the same in and out of the basal plane. 
Importantly, there is no evidence of any suppression of 
the gap in some direction or on some sheet of the Fermi 
surface. Indeed, as far as the quasiparticle transport is 
concerned, the superconducting gap appears to have the 
same uniform value everywhere on the Fermi surface. 

Relation to multi-band scenario. -The slow rise of kq/T 
at low H in LiFeAs is very different from the rapid rise 






LiFeAs 


H c2 (0) - 18 T 





10 12 

T C (K) 



14 16 18 20 



FIG. 3: Top panel: In-plane resistivity (p a ; small black dots) 
and Seebeck coefficient (thermopower) (S; large red dots) of 
LiFeAs, measured in zero magnetic field H, plotted as — S/T 
vs T. Both give a zero-field superconducting transition tem- 
perature T c — 18 K. The line is a quadratric fit to the p a (T) 
data below 50 K, extended to T = in order to extract 
an extrapolated value of the normal-state residual resistiv- 
ity po — 10 /x£l cm. The negative value of S indicates that 
electron-like carriers dominate the conductivity of LiFeAs at 
low temperature. Bottom panel: Temperature dependence of 
the superconducting upper critical field H C 2{T), determined 
by detecting T c in S/T vs T for different field strengths. The 
line is a smooth extrapolation to T = 0, giving an estimate 
of the zero-temperature critical field: H C 2(0) — 18 T. 



seen in typical multi-band superconductors such as MgB 2 
[39] and NbSe 2 [25] (see Fig. 2), in which the magni- 
tude of the s-wave superconducting gap is significantly 
different on two sheets of the Fermi surface. In both 
MgB 2 and NbSe 2 , the small gap is roughly one third 
of the large gap, which translates into the existence of 
a field scale H* ~ H C 2/9 sufficient to suppress super- 
conductivity on the small-gap Fermi surface, which can 
then contribute its full normal-state conductivity even 
deep inside the vortex state [231 I2S]- Specifically, at 
H = H c2 /5 > H*, k /T is already half (one third) of 
kn/T in MgB 2 (NbSe 2 ). If the gap on the electron Fermi 
surface of LiFeAs were 2 to 3 times larger than the gap 
on the hole Fermi surface, as reported by ARPES studies 
[12) . we would expect a significant enhancement of kq/T 
on a field scale H* ~ i? c2 /9 — _ff c2 /4. No such enhance- 
ment is observed. 

Two effects could possibly reconcile the small value of 
kq/T at low H in LiFeAs with a small gap on the hole 
Fermi surface. The first derives from the fact that it is 
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not the gap A that controls the tunneling, and hence 
the heat transport, but the coherence length fa °c vp/A 
[24] . A small value of vp on the hole surface could indeed 
compensate for the smaller gap. Specifically, if vp/v F = 
A e /Ah, then fa = fa and no multi-band feature in the 
H dependence of k /T is expected. ARPES data does 
suggests that vp > Vp [I2|T3] and it may be that fa ~ fa 
in LiFeAs. This would make H* ~ H c i- 

The second effect is if the normal-state conductivity of 
the hole Fermi surface were much smaller than that of the 
electron surface, i.e. if 07, << a e , or k'^/T << k f /T. 
The relative contribution of the small-gap hole Fermi sur- 
face at low H would then be a small fraction of the total 
(ko/T)/(ktsi/T) and hence difficult to resolve. Empirical 
evidence that < a e in LiFeAs comes from the fact 
that both Hall [10] and Seebeck (Fig. 3) coefficients are 
negative at low temperature. 

Comparison to other pnictides.- LiFeAs exhibits a 
temperature dependence of resistivity and a pressure de- 
pendence of T c that are consistent with an effective dop- 
ing level close to optimal (where T c is maximal) . Now, at 
optimal doping, both Co-Bal22 and K-Bal22 show a full 
isotropic gap in 3D 011011111112, just as reported here for 
LiFeAs. By contrast, in the low-T c stoichiometric super- 
conductors KFe2As2 pL[ and LaFePO [3HS] the supercon- 
ducting gap has nodes. This suggests that there may be a 
correlation between a high T c and a full, isotropic, node- 
less gap. In other words, high-temperature superconduc- 
tivity in iron-based materials would appear to thrive on 
an isotropic gap, in contrast with high-temperature su- 
perconductivity in copper oxides, which is intrinsically 
anisotropic and nodal. The only compound which shows 
nodal behavior at optimal doping is BaFe2(Asi_ x Pa;)2 
[7j. It remains to be seen whether this may be due to 
some unique feature of the multi-sheet Fermi surface in 
that material, such as a more pronounced c-axis disper- 
sion. 

Summary. - Our directional measurements of quasipar- 
ticle transport in the T = limit show that the supercon- 
ducting gap of LiFeAs is nodeless and isotropic in all di- 
rections. This excludes d-wave symmetry, and any other 
symmetry that requires line nodes on any piece of the 
multi-sheet Fermi surface of this superconductor. Sym- 
metries consistent with this constraint include s-wave 
and s± (whereby a full gap changes sign from the elec- 
tron Fermi surface to the hole Fermi surface [43]). A 
nodeless isotropic gap is also found in the iron-pnictidc 
superconductors Co-Bal22 [TO] and K-Bal22 [HI 1321 at 
optimal doping, suggesting a possible connection between 
isotropic gap and maximal T c . 
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